Test of a Jastrow-type wavefunction for a trapped few-body system in one dimension

For a system with interacting quantum mechanical particles in a one-dimensional harmonic oscillator, a trial wavefunction with simple structure based on the solution of the corresponding two-particle system is suggested and tested numerically. With the inclusion of a scaling parameter for the distance between particles, at least for the very small systems tested here the ansatz gives a very good estimate of the ground state energy, with the error being of the order of ~1% of the gap to the first excited state

Unrestricted Hartree-Fock theory of Wigner crystals

We demonstrate that unrestricted Hartree-Fock theory applied to electrons in a uniform potential has stable Wigner crystal solutions for $r_s \geq 1.44$ in two dimensions and $r_s \geq 4.5$ in three dimensions. The correlation energies of the Wigner crystal phases are considerably smaller than those of the fluid phases at the same density.Comment: 4 pages, 5 figure

Various spin-polarization states beyond the maximum-density droplet: a quantum Monte Carlo study

Using variational quantum Monte Carlo method, the effect of Landau-level mixing on the lowest-energy--state diagram of small quantum dots is studied in the magnetic field range where the density of magnetic flux quanta just exceeds the density of electrons. An accurate analytical many-body wave function is constructed for various angular momentum and spin states in the lowest Landau level, and Landau-level mixing is then introduced using a Jastrow factor. The effect of higher Landau levels is shown to be significant; the transition lines are shifted considerably towards higher values of magnetic field and certain lowest-energy states vanish altogether.Comment: 4 pages, 2 figures. Submitted to Phys. Rev.

Diffusion quantum Monte Carlo study of three-dimensional Wigner crystals

We report diffusion quantum Monte Carlo calculations of three-dimensional Wigner crystals in the density range r_s=100-150. We have tested different types of orbital for use in the approximate wave functions but none improve upon the simple Gaussian form. The Gaussian exponents are optimized by directly minimizing the diffusion quantum Monte Carlo energy. We have carefully investigated and sought to minimize the potential biases in our Monte Carlo results. We conclude that the uniform electron gas undergoes a transition from a ferromagnetic fluid to a body-centered-cubic Wigner crystal at r_s=106+/-1. The diffusion quantum Monte Carlo results are compared with those from Hartree-Fock and Hartree theory in order to understand the role played by exchange and correlation in Wigner crystals. We also study "floating" Wigner crystals and give results for their pair-correlation functions

Non-Equilibrium Electron Transport in Two-Dimensional Nano-Structures Modeled by Green's Functions and the Finite-Element Method

We use the effective-mass approximation and the density-functional theory with the local-density approximation for modeling two-dimensional nano-structures connected phase-coherently to two infinite leads. Using the non-equilibrium Green's function method the electron density and the current are calculated under a bias voltage. The problem of solving for the Green's functions numerically is formulated using the finite-element method (FEM). The Green's functions have non-reflecting open boundary conditions to take care of the infinite size of the system. We show how these boundary conditions are formulated in the FEM. The scheme is tested by calculating transmission probabilities for simple model potentials. The potential of the scheme is demonstrated by determining non-linear current-voltage behaviors of resonant tunneling structures.Comment: 13 pages,15 figure

Optimization of inhomogeneous electron correlation factors in periodic solids

A method is presented for the optimization of one-body and inhomogeneous two-body terms in correlated electronic wave functions of Jastrow-Slater type. The most general form of inhomogeneous correlation term which is compatible with crystal symmetry is used and the energy is minimized with respect to all parameters using a rapidly convergent iterative approach, based on Monte Carlo sampling of the energy and fitting energy fluctuations. The energy minimization is performed exactly within statistical sampling error for the energy derivatives and the resulting one- and two-body terms of the wave function are found to be well-determined. The largest calculations performed require the optimization of over 3000 parameters. The inhomogeneous two-electron correlation terms are calculated for diamond and rhombohedral graphite. The optimal terms in diamond are found to be approximately homogeneous and isotropic over all ranges of electron separation, but exhibit some inhomogeneity at short- and intermediate-range, whereas those in graphite are found to be homogeneous at short-range, but inhomogeneous and anisotropic at intermediate- and long-range electron separation.Comment: 23 pages, 15 figures, 1 table, REVTeX4, submitted to PR

Metal Surface Energy: Persistent Cancellation of Short-Range Correlation Effects beyond the Random-Phase Approximation

The role that non-local short-range correlation plays at metal surfaces is investigated by analyzing the correlation surface energy into contributions from dynamical density fluctuations of various two-dimensional wave vectors. Although short-range correlation is known to yield considerable correction to the ground-state energy of both uniform and non-uniform systems, short-range correlation effects on intermediate and short-wavelength contributions to the surface formation energy are found to compensate one another. As a result, our calculated surface energies, which are based on a non-local exchange-correlation kernel that provides accurate total energies of a uniform electron gas, are found to be very close to those obtained in the random-phase approximation and support the conclusion that the error introduced by the local-density approximation is small.Comment: 5 pages, 1 figure, to appear in Phys. Rev.

Electronic structure of rectangular quantum dots

We study the ground state properties of rectangular quantum dots by using the spin-density-functional theory and quantum Monte Carlo methods. The dot geometry is determined by an infinite hard-wall potential to enable comparison to manufactured, rectangular-shaped quantum dots. We show that the electronic structure is very sensitive to the deformation, and at realistic sizes the non-interacting picture determines the general behavior. However, close to the degenerate points where Hund's rule applies, we find spin-density-wave-like solutions bracketing the partially polarized states. In the quasi-one-dimensional limit we find permanent charge-density waves, and at a sufficiently large deformation or low density, there are strongly localized stable states with a broken spin-symmetry.Comment: 8 pages, 9 figures, submitted to PR

An Effective-Medium Tight-Binding Model for Silicon

A new method for calculating the total energy of Si systems is presented. The method is based on the effective-medium theory concept of a reference system. Instead of calculating the energy of an atom in the system of interest a reference system is introduced where the local surroundings are similar. The energy of the reference system can be calculated selfconsistently once and for all while the energy difference to the reference system can be obtained approximately. We propose to calculate it using the tight-binding LMTO scheme with the Atomic-Sphere Approximation(ASA) for the potential, and by using the ASA with charge-conserving spheres we are able to treat open system without introducing empty spheres. All steps in the calculational method is {\em ab initio} in the sense that all quantities entering are calculated from first principles without any fitting to experiment. A complete and detailed description of the method is given together with test calculations of the energies of phonons, elastic constants, different structures, surfaces and surface reconstructions. We compare the results to calculations using an empirical tight-binding scheme.Comment: 26 pages (11 uuencoded Postscript figures appended), LaTeX, CAMP-090594-

Computational Physics on Graphics Processing Units

The use of graphics processing units for scientific computations is an emerging strategy that can significantly speed up various different algorithms. In this review, we discuss advances made in the field of computational physics, focusing on classical molecular dynamics, and on quantum simulations for electronic structure calculations using the density functional theory, wave function techniques, and quantum field theory.Comment: Proceedings of the 11th International Conference, PARA 2012, Helsinki, Finland, June 10-13, 201